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ABSTRACT 

We present the supernova (SN) sample and Type-la SN (SN la) rates from the Cluster Lensing And Supernova 
survey with Hubble (CLASH). Using the Advanced Camera for Surveys and the Wide Field Camera 3 on the 
Hubble Space Telescope (HST), we have imaged 25 galaxy-cluster fields and parallel fields of non-cluster galaxies. 
We report a sample of 27 SNe discovered in the parallel fields. Of these SNe, ~ 13 are classified as SN la candidates, 
including four SN la candidates at redshifts z > 1.2. We measure volumetric SN la rates to redshift 1.8 and add the 
first upper limit on the SN la rate in the range 1 .8 < z < 2.4. The results are consistent with the rates measured by 
the HST/ GOODS and Subaru Deep Field SN surveys. We model these results together with previous measurements 
at z < 1 from the literature. The best-fitting SN la delay-time distribution (DTD; the distribution of times that 
elapse between a short burst of star formation and subsequent SN la explosions) is a power law with an index of 
- l.OO^S 09 ', (statistical)! 0 ; ^ (systematic), where the statistical uncertainty is a result of the 68% and 95% (in 
parentheses) statistical uncertainties reported for the various SN la rates (from this work and from the literature), 
and the systematic uncertainty reflects the range of possible cosmic star-formation histories. We also test DTD 
models produced by an assortment of published binary population synthesis (BPS) simulations. The shapes of 
all BPS double-degenerate DTDs are consistent with the volumetric SN la measurements, when the DTD models 
are scaled up by factors of 3-9. In contrast, all BPS single-degenerate DTDs are ruled out by the measurements 
at > 99% significance level. 

Key words: supernovae: general - surveys - white dwarfs 
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E INTRODUCTION 

Although Type-la supernovae (SNe la) have been used to 
measure extragalactic distances and thus reveal the accelerating 
expansion of the universe (Riess et al. 1998; Schmidt et al. 
1998; Perlmutter et al. 1999), the nature of the stellar system 
that leads to these explosions remains unclear (see review by 
Howell 201 1). The current consensus is that the progenitor is a 
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carbon-oxygen white dwarf (CO WD) that accretes matter from 
a binary companion until the pressure or temperature somewhere 
in the WD become high enough to ignite the carbon and lead 
to a thermonuclear explosion of the WD (Leibundgut 2000). 
Different scenarios have been proposed to explain the nature 
of the binary companion and the process of mass accretion. 
The leading scenarios are the single-degenerate (SD) scenario 
(Whelan & Iben 1 973), in which the binary companion is either a 
main-sequence star, a subgiant just leaving the main sequence, a 
red giant, or a stripped “He star,” and the WD accretes mass from 
the secondary through Roche-lobe overflow or a stellar wind. 


1 


The Astrophysical Journal , 783:28 (19pp), 2014 March 1 


Graur et al . 


In the double-degenerate (DD) scenario (Iben & Tutukov 1984; 
Webbink 1984), the companion is a second CO WD and the two 
WDs merge due to loss of energy and angular momentum to 
gravitational waves. 

Each of these scenarios predicts a different fonn of the 
distribution of times that elapse between a short burst of star 
formation and any subsequent SN la events, known as the delay- 
time distribution (DTD; see Wang & Han 2012 and Hillebrandt 
et al. 2013 for recent reviews). The DTD can be thought of as 
a transfer function connecting the star- formation history (SFH) 
of a specific stellar environment and that environment’s SN la 
rate. Thus, by measuring the SN la rate and comparing it to 
the SFH, one might reconstruct the DTD. The SN la DTD has 
been recovered using several techniques applied to different SN 
samples collected from different types of stellar environments 
(see review by Maoz & Mannucci 2012). The emerging picture 

is that of a power-law DTD with an index of 1, a form that 

arises naturally from the DD scenario, although combinations 
of DTDs from a DD channel and a SD channel cannot be ruled 
out. One method to recover the DTD, Tkt), is to measure the SN 
la rate, Ri a (t), as a function of cosmic time t in field galaxies, 
and compare them to the cosmic SFH, S(t): 

R Ia (t)= S(t-xWT)dT. (1) 

0 

Measurements of the volumetric SN la rates (i.e., the SN la 
rates per unit volume) in field galaxies agree out to z = 1. 
Graur et al. (2011, G1 1) provide a compilation of all SN la rates 
measured up to 201 1, and later measurements are presented by 
Krughoffet al. (201 1), Perrett et al. (2012), Barbary et al. (2012), 
Melinder et al. (2012), and Graur & Maoz (2013). Volumetric 
SN la rate measurements were first extended to z > 1 by Dahlen 
et al. (2004), with additional data analyzed by Dahlen et al. 
(2008, D08), using the Hubble Space Telescope (HST) to survey 
the GOODS fields (Giavalisco et al. 2004; Riess et al. 2004). 
The HST/ GOODS survey discovered 20 SNe la at 1 < z < 1.4 
and 3 at 1.4 < z < 1.8. Gil conducted a SN survey in the 
Subaru Deep Field (SDF) using the 8.2 m Subaru telescope 
and discovered 27 SN la candidates at 1 < z < 1.5 and 10 at 
1.5 < z < 2. 

The SN la rate uncertainties at z > 1, and especially at 
z > 1.5, are dominated by small-number statistics. The three 
z > 1.4 HST/ GOODS SNe la were discovered in host galaxies 
having a spectroscopic redshift (spec-z) and no active galactic 
nucleus (AGN) activity. On the other hand, the classification 
of the larger SDF SN sample from Gil relies on photometric 
redshift (photo-z) measurements that might be systematically 
biased toward high redshifts (see their Section 4.2). While Gil 
also used several methods to weed out interloping AGNs, there 
could still be some AGN contamination because each SN in 
the SDF sample was only observed on one epoch. The HST/ 
GOODS sample, while smaller than the SDF sample, suffers 
from lower systematic uncertainties owing to the spectroscopic 
classification of the SN host galaxies and measurements of their 
redshifts, and a better sampling of the SN light curves. Of the 10 
z > 1.5 SN host galaxies in the SDF sample, only one galaxy 
has so far had its redshift and lack of AGN activity confirmed 
spectroscopically (Frederiksen et al. 2012). 

Although the GOODS and SDF z > 1 SN la rates are 
consistent, their interpretation differs between the two groups. 
Based on the GOODS data, Dahlen et al. (2004, 2008) argued 
that the SN la rate declined at z > 0.8. Fitting this declining SN 
la rate evolution, Strolger et al. (2004, 2010) surmised that the 


DTD is confined to delay times of 3-4 Gyr. In contrast, based 
on the SDF data, Gil found that the SN la rate evolution does 
not decline at high redshifts, but rather levels off, as would be 
expected of a power-law DTD. 

Two new SN surveys are attempting to resolve this conflict. 
These surveys are components of two 3 yr HST Multi-Cycle 
Treasury programs that use the Advanced Camera for Surveys 
(ACS) and the new Wide Field Camera 3 (WFC3). Results 
from the Cosmic Assembly Near-infrared Deep Extragalactic 
Legacy Survey (CANDELS; Grogin et al. 2011; Koekemoer 
et al. 2011) will be reported by Rodney et al. (2014). Here, 
we describe results from the Cluster Lensing And Supernova 
survey with Hubble (CLASH; Postman et al. 2012). CLASH 
imaged 25 galaxy clusters in 1 6 broad-band filters from the near- 
ultraviolet (NUV) to the near-infrared (NIR) with the ACS and 
WFC3 cameras working in parallel mode. While one camera 
was pointed at the galaxy cluster, the other one was used to 
observe a parallel field far enough from the galaxy cluster so as 
not to be significantly affected by strong lensing. 

In this work, we report a sample of 27 SNe discovered 
in the parallel fields of the 25 CLASH galaxy clusters. In 
Section 2, we describe the CLASH observations and our 
imaging and spectroscopic follow-up program. We report our 
SN sample in Section 3, where we also conduct detection- 
efficiency simulations and classify the SNe. Using our SN la 
sample, we measure SN la rates out to z “ 2.4 in Section 4 
and use them to test different forms of the DTD in Section 5. 
Finally, we summarize our results in Section 6. Throughout 
this work, we assume a A-cold-dark-matter cosmological model 
with parameters Qa = 0.7, Q m = 0.3, and H 0 = 70 km s" 1 
Mpc" 1 . Unless noted otherwise, all magnitudes are on the Vega 
system. 

We designate our SN candidates according to the cluster and 
year in which they were discovered and the first three letters of 
the nickname given to them for internal tracking purposes. For 
example, CLI1 lHad is a CLASH (CL) SN that was discovered 
in one of the parallel fields around the ninth (or Ith) cluster, 
MACS0717. 5+3745, in 2011, and was nicknamed “Hadrian.” 
For the sake of brevity, we will henceforth refer to our SN 
candidates simply as SNe. 

2. OBSERVATIONS 
2.1. Imaging 

The CLASH observation strategy is described in detail by 
Postman et al. (2012). During Cycles 18-20, CLASH observed 
25 galaxy clusters in the redshift range 0.187-0.890. The central 
region of each galaxy cluster was imaged with 16 broad-band 
filters from the NUV to the NIR using the ACS and WFC3 
cameras on HST. In ACS, we used the Wide Field Channel 
(WFC), with a field of view of 202 x 202 and a pixel scale of 
0. 05 pixel" 1 . WFC3 includes two detectors: an infrared channel 
(WFC3-IR) with a field of view of 123 * 136 and a scale of 
0. 13 pixel" *; and an ultraviolet-visible channel (WFC3-UVIS) 
with a field of view of 162 * 162 and a scale of 0. 04 pixel" 1 . 

The orientation of HST and the cadence between succeed- 
ing visits to the galaxy cluster (“prime”) field were chosen so 
that two ACS and two WFC3 parallel fields would each be ob- 
served on four separate occasions, with a median cadence of 
18 days. Each visit to a WFC3 parallel field consisted of one 
orbit comprising two F160W filter exposures and one exposure 
in filters F125W and F350LP each (filter+system central wave- 
lengths X 0 ~ 15,369, 12,486, and 5846 A, respectively). Visits 
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Table 1 

Typical Exposure Times for CLASH Parallel Fields 


Camera 

Filter 

Exposures 

Total Time 
(s) 

5a Limiting Magnitude 
(Vega mag) 

ACS-WFC 

F850LP 

4 

1500 

25.0 



6 

3600 

25.4 

ACS-WFC 

F775W 

1 

400 

25.7 



2 

700 

25.9 

WFC3-IR 

F160W 

2 

1200 

25.4 

WFC3-IR 

F125W 

1 

700 

25.7 

WFC3-UVIS 

F350LP 

1 

650 

27.5 


to the ACS parallel fields consisted of one orbit when the prime 
field was imaged with either the ACS or WFC3-IR cameras and 
two orbits when the prime field was imaged with the WFC3- 
UV1S camera. During single-orbit visits, the parallel ACS orbit 
comprised four F850LP filter exposures and one F775W fil- 
ter exposure (filter+system central wavelengths /, (l = 9445 and 
7764 A, respectively). When the ACS parallel fields were im- 
aged over two orbits, they consisted of six F850LP and two 
F775W exposures. These filters, the reddest in each camera, 
were chosen to detect high-redshift SNe. The F350LP band was 
added to the WFC3 observations for additional color informa- 
tion to aid in the classification of any SNe discovered in those 
fields. The HST angular resolution in our search bands is ~ 0. 10 
and —0.17 in F850LP and F160W, respectively, slightly larger 
than the pixel scales of their respective cameras. Table 1 lists 
the typical exposure times and 5a limiting magnitudes reached 
in each of these filters. 

The limiting magnitude in each filter was calculated using 
the method outlined in Kashikawa et al. (2004): we conducted 
aperture photometry on hundreds of blank regions in the image, 
fit a Gaussian to the negative side of the resultant histogram (as 
the positive tail could be contaminated by light from the sources 
in the image), and treated the standard deviation of the fit as 
an estimate of the average noise in the image. We used circular 
apertures with radii of 4, 3, and 5 pixels, which correspond 
to 0. 20, 0. 27, and 0. 20 in the pixel scale to which we drizzle 
the images taken with ACS, WFC3-IR, and WFC3-UVIS: 0.05, 
0.09, and 0.04 arcsec pixel - ', respectively. 

An additional 52 HST orbits were allocated for follow- 
up imaging or slitless spectroscopy of targets of opportunity, 
such as high-redshift SNe la. This cache of orbits was added 
to the 150 similar HST orbits allocated to the CANDELS 
program, for a sum of 202 follow-up orbits for the combined 
CLASH+CANDELS SN survey (PI: A. Riess). 

Our HST reduction and image-subtraction pipeline is de- 
scribed in detail by Rodney et al. (2014). Briefly, the raw HST 
images were first calibrated using the STSDAS 22 calibration 
tools. The calibrations include bias correction, dark subtraction, 
and flat fielding. In the case ofWFC3-lR images, “up-the-ramp” 
fitting was used to remove cosmic ray (CR) events. Charge- 
transfer efficiency losses in the ACS images were corrected us- 
ing the algorithm of Anderson & Bedin (2010). Next, the subex- 
posures in each filter were combined using Mul t iDrizzl e 
(Koekemoer et al. 2003). This stage also removed the geometri- 
cal distortion of the HST focal plane. For each filter, we created 
“template” images comprised of all previous observations in the 
same filter. This means that some SN light may be included in 
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the template images, which we take into account in Section 4. 
Finally, we subtracted the template images from the drizzled 
“target” images to produce the difference images that were then 
searched for SNe. Owing to the stable point-spread-function 
(PSF) of HST, we did not need to degrade the PSF of either the 
target or template images to match the PSF of the images, as 
done in ground-based SN surveys (e.g., Gil). 

Most of the CLASH galaxy clusters were observed in the B, V, 
Rc, I c , and z bands with Suprime-Cam (Miyazaki et al. 2002) at 
the prime focus of the 8.2 m Subaru telescope, for the purpose of 
measuring the amount of shear induced on background galaxies 
by weak lensing from the galaxy cluster, and for deriving the 
photometric redshifts of the galaxies in the CLASH parallel 
fields. For an example of such observations of the CLASH 
galaxy cluster MACS 1206, and a description of their reduction, 
see Umetsu et al. (2012). 

2.2. Spectroscopy 

The host galaxies of all SN candidates, and in several cases 
the SNe themselves, were followed up with spectroscopic ob- 
servations from several ground-based observatories, as detailed 
below, or with HST slitless spectroscopy, using the G800L ACS 
grism spectrograph. The ground-based observatories and in- 
struments used for this work were the Low Resolution Imaging 
Spectrometer (LRIS; Oke et al. 1995) and the DEep Imaging 
Multi-Object Spectrograph (DEIMOS; Faber et al. 2003) on 
the Keck 1 and II 10 m telescopes, respectively; the Gemini 
Multi-Object Spectrograph (GMOS; Hook et al. 2003) on the 
Gemini North and South telescopes (GeminiN and GeminiS, 
respectively); the Multi-Object Double Spectrograph (MODS; 
Pogge et al. 2010) on the Large Binocular Telescope (LBT); and 
the FOcal Reducer and low dispersion Spectrograph (FORS; 
Appenzeller et al. 1998), the Visible MultiObject Spectrograph 
(VIMOS; Le Fevre et al. 2003), and the X-shooter spectro- 
graph (Vernet et al. 2011) on the Very Large Telescope (VLT). 
Table 2 details which instruments were used to obtain spectra of 
each SN host galaxy. Several examples of SN host-galaxy spec- 
tra are shown in Figure 1. The spectra of the SNe CLFllVes, 
CLIllHad, and CLY13Pup are shown in Figure 2. 

3. SUPERNOVA SAMPLE 

In our survey, SNe can be discovered if they either brighten 
or decline between one search epoch and the next. By a 
“brightening” SN we do not mean that the SN is necessarily 
caught on the rising part of its light curve, but rather any case in 
which the discovery flux is higher than the template flux. As a 
result of the cadence of our survey, it is easier to discover SNe 
either on the rise or near peak, as in these cases the template 
image will contain no SN flux. In contrast, SNe caught while 
on the decline will invariably have some flux in all our images, 
thus reducing the flux in the difference image and consequently 
their probability of detection. We have discovered a total of 20 
brightening and 7 declining SNe in the parallel fields of the 25 
CLASH clusters. Of these, 18 were discovered in the ACS and 
9 in the WFC3 fields. Nineteen (or 70%) of the SN host galaxies 
have spectroscopic redshifts, as detailed below in Section 3.3. 
We classify half of this sample as SNe la, four of which are at 
z > 1.2. Our SN sample is summarized in Table 2 and displayed 
in F igure 3 . 

We have discovered 12 additional SN candidates in the prime 
fields. However, as the effects of gravitational lensing must 
be taken into account to properly classify any SNe discovered 
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Table 2 

SNe Discovered in the Parallel Fields of the 25 CLASH Galaxy Clusters 


ID 

Nickname 

a (h m s) 

5f ) 

P (Ia)wp 

P (I a )np 

Photo-z 

Spec-z 

Spec-z Source 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

CLAlOCal 

Caligula 

02 : 48 : 25.740 

-03 : 33 : 08.37 

Q qc+ 0.03 
U ‘ y:> -0.14 

0-98-q.oo 

1-68^0 *15 


(VLT+X-shooter) 

CLAlONer 

Nero 

02 : 47 : 40.180 

-03 : 32 : 53.29 

0 76 +0 09 
u - ' - 0.26 

0- 82^0.08 

O.32! 0 o°o 8 1 

0.362 

Keck+DEIMOS 

CLBllOth 

Otho 

1 1 : 49 : 56.745 

+22 : 18 : 42.87 

0 70+0.06 
u - /y -0.14 

0 89 +0 00 
u - oy - 0.01 

0 - 86 ^ 0.20 

0.962 

VLT+FORS2 

CLCllTit 

Titus 

17 : 22 : 52.985 

+32 : 07 : 25.74 

0 89 +003 
u - oy -0.08 

0 94 + 0.01 
u 0.01 

0 70 +0 - 05 
u - ' -0.05 

0.839 

Keck+LRIS 

CLDllCla 

Claudius 

12 : 06 : 08.868 

-08 : 42 : 54.73 

0 i -1+0.24 
U - 1J -0.13 

0 19 +O.I 8 

y .lV_o 19 

0 24 +0 - 07 
UlZ ^-0.04 


(VLT+FORS2; Keck+LRIS; ACS+G800L) 

CLEllAug 

Augustus 

13 : 47 : 12.802 

- 11 : 42 : 28.97 

0 OJ +0 06 
U-U ' ) -0.03 

0 06 +0 05 
u - uo -0.06 

0 714 + 0.10 
u,:> ^-0.03 

0.329 

GeminiS+GMOS 

CLF llVes 

Vespasian 

21 : 29 : 42.612 

-07 : 41 : 48.08 

0 97 + 0.01 
u - y/ -0.04 

0 98 +0 00 
u - y °-o.oi 

111 +0-05 

1 - 31 -o.io 

1.22 

ACS+G800L; Keck+DEIMOS 

CLFllDom 

Domitian 

21 : 29 : 53.224 

-07 : 40 : 56.95 

n £1+0. 16 
U.o3_0 40 

0 7 1 + 0.06 
U. 7 1-0.24 

0 71 + 0.14 

u - ^-o.oo 



CLHllTra 

Trajan 

21 : 39 : 46.036 

-23 : 38 : 34.71 

i oo +0 0 ° 
A>uu - 0.00 

i oo +0 0 ° 
AUU - 0.00 

1 41 +014 
A -^ A - 0.11 


(VLT+VIMOS; Keck+DEIMOS) 

CLlllHad 

Hadrian 

07 : 17 : 20.115 

+37 : 49 : 54.54 

i oo +0 0 ° 
A - uu -o.oo 

i oo +0 0 ° 
AUU - 0.00 

0 2 1 +0 - 03 
u - zl -0.03 

0.261 

Keck+LRIS 

CLKllBur 

Burgundy 

06 : 49 : 13.878 

+70 : 13 : 17.91 

o oo +00 ° 
u - uu -o.oo 

o oo +00 ° 
u - uu -o.oo 

0 27 +0 05 
- 0.05 

0.281 

Keck+LRIS 

CLIllPiu 

Antoninus Pius 

07 : 17 : 59.057 

+37 : 40 : 51.13 

o oo +00 ° 
u - uu -o.oo 

o oo +00 ° 
uuu - 0.00 

A 1 O+ 0.10 
U - A °- 0.02 

0.191 

GeminiS+GMOS; LBT+MODS 

CLL12Aur 

Marcus Aurelius 

11 : 16 : 07.758 

+01 : 23 : 44.60 

0 00 + 0.00 
u - uu -o.oo 

0 00 + 0.00 
u - uu -o.oo 

0 79 + 0.10 
u,zy -0.04 

0.271 

Keck+DEIMOS 

CLL12Luc 

Lucius Verus 

11 : 15 : 57.141 

+01 : 23 : 19.62 

0 00 + 0.00 
u - uu -o.oo 

0 00 + 0.00 
u - uu -o.oo 

A 77 +0.06 

0.03 

0.36 

Keck+DEIMOS 

CLM12Com 

Commodus 

08 : 00 : 52.385 

+36 : 07 : 37.31 

o oo +00 ° 

u - uu -o.oo 

o oo +00 ° 
uuu - 0.00 

0 19 + 0.05 
u - iy -0.07 

0.207 

Keck+DEIMOS 

CLM12Car 

Cardinal 

08 : 00 : 42.802 

+36 : 07 : 13.81 

0 22 +0 - 27 
w - zz - 0.22 

0 54 +O.I 6 
0.34 

0 47 + 0.08 
u - 4/ -0.05 

0.518 

Keck+DEIMOS 

CLP12Get 

Geta 

21 : 29 : 23.918 

+00 : 08 : 24.77 

i oo +00 ° 
1 - uu -o.oo 

i oo +0 0 ° 
A - uu -o.oo 

1 64 +004 
a -°^-0.04 

1.64 

VLT+X-shooter; (VLT+VIMOS) 

CLR12Arm 

Neill Armstrong 

01 : 31 : 30.231 

- 13 : 34 : 39.13 

0 0 1 + 0.02 
u - ul -o.oi 

0 01 + 0.01 
u - UA -o.oi 

1 1 7+0.63 
a - 1z -0.36 


(VLT+VIMOS; Keck+LRIS) 

CLS12Mac 

Macrinus 

04 : 15 : 47.671 

-24 : 00 : 23.77 

0 O9 +016 
u - uy -0.09 

0 75 +O. 13 
u - ZJ - 0.22 

0 95 + 0.08 
u - y3 -0.03 

1.034 

Keck+DEIMOS 

CLT 1 2Ela 

Elagabalus 

22 : 48 : 09.132 

-44 : 35 : 16.63 

0 01+0.06 
u - ol - 0.10 

0 70 +0 07 
u - ' - 0.06 

0 C 5+0.03 
u - J - ;, -0.04 

0.6058 

Keck+LRIS 

CLT 1 2Ale 

Alexander Severus 

22 : 49 : 20.961 

-44 : 32 : 47.94 

0 12 +018 
u - AZ -o.n 

o n +01 ° 
u - 1 :) -o.io 

1 o +2 ° 

0.5 


(VLT+FORS2) 

CLC12Thr 

Thrax 

17 : 22 : 44.529 

+32 : 03 : 35.96 

o oo +00 ° 
u - uu -o.oo 

o oo +00 ° 
uuu - 0.00 

0 2 J +0 - 03 
u,ZJ -0.04 



CLV12Gor 

Gordian 

11 : 57 : 10.258 

+33 : 42 : 19.60 

0 00 + 0.00 
u - uu -o.oo 

0 00 + 0.00 
uuu - 0.00 

0 57 +O.I 8 

U - J 7 — 0. 14 

0.516 

Keck+DEIMOS 

CLY13Lil 

Lilia 

13 : 11 : 03.822 

-03 : 15 : 50.66 

i oo +0 0 ° 

A - UL, -o.oo 
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Notes. (1) SN identification. (2) SN nickname. (3)-(4) Right ascensions and declinations (J2000). (5)-(6) STARDUST probability of classification as a SN la with 
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behind the galaxy clusters, we leave their treatment to a future 
paper. The complete photometry of all 39 SNe in our sample will 
appear in a future paper by O. Graur et al. (2014, in preparation). 

3.1. Candidate Selection 

The F850LP- and F160W-band subtraction images were 
simultaneously searched by eye and scanned with the source- 
identifying software SEXTRACTOR (Bertin & Arnouts 1996) 
to identify variable objects. We set SEXTRACTOR to locate all 
objects that had at least four connected pixels with flux 3o above 
the local background level in both the regular subtraction image 
and in its negative, the latter in order to search for declining 
SNe. To increase the detection efficiency in the F160W band, 
the FI 60 W- andF125W-band subtraction images were searched 
by eye simultaneously (by toggling between them), as some SNe 
may appear brighter in the F125W band (see, for example, the 
light-curve fit of CLlllHad in Figure 6). The F850LP- and 
F775W-band subtraction images were not similarly toggled due 
to the high CR contamination in the F775W-band subtraction 
image. 

To be regarded as SN candidates, the variable objects had to 
pass the following criteria. 


1. All objects with suspect residual shapes, such as the 
subtraction residuals of bright galaxy cores or objects with 
non-PSF shapes, were rejected. 

2. The F850LP- and F160W-band images were comprised 
of several subexposures (four or six in the F850LP band 
and two in the F160W band). We used these subexposures 
to create separate subtraction images, which were then 
compared to the main subtraction image. The object had 
to appear in all of the subtraction images to be considered 
a likely candidate. 

3. To be considered a declining SN, the object had to have a 
negative flux in the subtraction image and appear in both the 
target and template images. Objects that only appeared in 
the template images were discarded as either CRs or noise 
spikes. 

4. Objects with suspect residual shapes in the F775W- or 
F125W-band subtraction images were flagged for inspec- 
tion in the next search epoch. As a result of the cadence of 
the survey, any SNe detected in one of the first two search 
epochs would be visible in the other search epochs as well, 
so if the object under consideration did not reappear in a 
later epoch, it was rejected. No candidates were rejected if 
they did not appear in the F775W- or F125W-band images. 
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Figure 1. Examples of SN host-galaxy spectra. From top to bottom, we present 
spectra of the SN host galaxies CLI1 lHad, atz = 0.261, taken with Keck+LRIS; 
CLEllAug, at z = 0.329, taken with GeminiS+GMOS (the gaps in this 
spectrum are the result of physical gaps between GMOS chips); CLCllTit, at 
z = 0.839, obtained with Keck+LRIS; CLS12Mac, at z = 1.034, obtained with 
Keck+DEIMOS; and CLFllVes, at z = 1.22, obtained with Keck+DEIMOS. 
All spectra have been binned into 10 A wide bins. All flux units and have been 
arbitrarily scaled. 

3.2. Detection-efficiency Simulations 

In our survey, SNe can be missed because of many factors. 
Generally, the fainter the SN, the less likely it is to be detected 
above the background. On average, F775W-band images suffer 
from a background (composed of zodiacal light, earthshine, 
and airglow) level twice as high as F850LP-band images 
(Sokol et al. 2012). The main sources of background for 
WFC3-IR observations are earthshine and zodiacal light, with 
the latter being the dominant source. Both sources contribute 
less background at longer wavelengths (Giavalisco et al. 2002), 
but as our F160W exposures are roughly twice as long as 
the F125W exposures, they both display roughly the same 
number of background counts (Dressel 2012). There are other 
factors that affect the discovery probability of a SN, such as its 
proximity to the core of its host galaxy (SNe that explode close 
to the cores of their host galaxies are harder to discover due to 
the noise from the higher background and the residuals from 
imperfect image subtractions). 

To test the effect of these and other factors on our detection 
efficiency, we planted ~ 1000 fake point sources in the raw 
images at the start of our reduction pipeline. The fake SNe 
were planted in random locations around galaxies chosen from 
SEXTRACTOR catalogs of the images following a Gaussian 
distribution centered on the center of the galaxy, as measured by 
SEXTRACTOR, with a standard deviation of o = 2 R50, where 



3000 4000 5000 6000 7000 8000 9000 10000 

Observed wavelength [A] 


Figure 2. SN spectra. In black solid curves, we show the spectra of the SNe 
CLIllHad (top), at z = 0.261, taken with Keck+LRIS; CLY13Pup (center), 
at z = 0.804, taken with VLT+FORS2; and CLFllVes (bottom), at z = 1.22, 
taken with the HST G800L ACS grism spectrograph. Overlaid on the spectra as 
red dashed curves are examples of SNe la from the literature shifted to the same 
redshift and scaled so as to fit the data. The SN designation, along with that of 
the example SN la, are noted on the top left of each panel. The top and center 
spectra have been binned into 1 0 A wide bins, while the bottom spectrum has 
been binned into 80 A wide bins. All flux units are arbitrary. Beyond 5500 A, the 
spectrum of CLI 1 1 Had may be dominated by host-galaxy light, as no correction 
for host-galaxy light was performed during reductions. 

(A color version of this figure is available in the online journal.) 

R50 is the radius that contains 50% of the galaxy light. This 
distribution assured that the fake SNe approximately followed 
the light of the galaxy and that a large number were planted 
in galaxy cores (e.g., Forster & Schawinski 2008). Near the 
center of a bright galaxy, a SN could be obscured due to the 
increased Poisson noise and residual subtraction artifacts from 
small inter-epoch registration errors. This is evaluated in Rodney 
et al. (2014), where we find this effect to be negligible, with less 
than 2% of galaxies above z = 0.2 exhibiting core residuals 
that could obscure a SN. The magnitudes of the fake SNe were 
drawn from flat distributions inF850LP andF160W in the range 
22-28 mag. To simulate the appearance of real SNe la, the 
F775W and F125W magnitudes, respectively, were randomly 
chosen from a SN la simulation done with the SuperNova 
ANAlysis (SNANA 23 ; Kessler et al. 2009b) software package, 
which was constructed to reflect a realistic spread of SN la 
colors in the redshift range z = 0-3, with host-galaxy extinction 
according to values chosen from an exponential of the form 
P (A v ) = e ( ~ Av/ ' tv) , with t v = 0.7, chosen to approximate the 
host-galaxy extinction model of Riello & Patat (2005). This was 
done to ensure that the fake SNe resembled the colors of their 
real counterparts as close as possible, and was of importance 
mainly in the WFC3 fields, where the SN searchers toggled 
between the F160W and F125W difference images. The PSF 
was simulated usingTiny Tim 24 (Krist et al. 2011). 

Figure 4 shows our detection efficiency, as a function of 
the brightness of the fake SNe, in the F850LP and F160W 
bands. The uncertainties of the measurements represent the 68% 
binomial confidence intervals. We follow Sharon et al. (2010) 


23 http://sdssdp62.fnal.gov/sdsssn/SNANA-PUBLIC/ 

24 http://www.stsci.edu/hst/observatory/focus/TinyTim 
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Figure 3. SNe discovered in the parallel fields of the CLASH clusters. North is up and east is left. In the triplet of tiles for each event, the left-hand tiles show the SN 
host galaxies without any SN light, whereas the center tiles display the SN host galaxy as imaged when the SN was first discovered. For the declining SNe CLK1 lBur, 
CLL12Luc, CLAlONer, CLV12Gor, CLF1 lDom, CLT12Ela, and CLY13Gal, the left-hand and center tiles show the SN and host galaxy on the first and last visits to the 
field, respectively. The right-hand tiles show the subtraction in the F850LP or F160W bands for SNe discovered in the ACS or WFC3 parallel fields, respectively. The 
stretch of the images and the location of the SN differ from panel to panel in order to highlight host-galaxy properties. The header of each panel gives the designation 
of the SN along with its redshift and camera. Spectroscopic redshifts (cases with no uncertainties in z noted) are given to three significant digits. Photometric redshifts 
are shown with their uncertainty; in cases where the photometric redshift is not well constrained, we note the approximate peak of the probability density function. 
(A color version of this figure is available in the online journal.) 
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Figure 4. SN detection efficiency vs. magnitude in the F850LP (left) and FI 60 W (right) bands. The uncertainties of the measurements are the 68% binomial confidence 
intervals. The dotted lines mark where the best-fit efficiency curves drop to 50%, at 25.2 and 25.0 mag in the F850LP and F160W bands, respectively. 


and fit the efficiency measurements with the function 

m-m l/2 ^ 

1 + e s i , m mi /2 

n(m;m 1/2 ,s 1 ,s 2 )= , m _ mi/2 -i (2) 

i 1 + e s 2 , m> mi/ 2 , 

where m is the magnitude in the F850LP band; mi /2 is the 
magnitude at which the efficiency drops to 50%; and Si and 
S 2 determine the range over which the efficiency drops from 
100% to 50%, and from 50% to 0, respectively. Our detection 
efficiency drops to 50% at 25.2 and 25.0 mag in the F850LP and 
F160W bands, respectively. 

3.3. Host-galaxy Redshifts 

Our classification method, as with most SN classification 
techniques, relies on a good knowledge of the redshift of either 
the SN or its host galaxy. As part of our survey strategy, we 
have endeavored to obtain spec-z measurements of the host 
galaxies of all the SNe in our sample, mostly with ground-based 
observatories, as described above. Some of the SN host galaxies 
suspected of being at z > 1.2 (CLD1 lCla and CLF1 IVes) were 
also followed up with HST slitless spectroscopy using the ACS 
G800L grism. At this time, we have acquired and reduced the 
spectra of 19 of the 27 SN host galaxies in our sample. 

For the remaining eight SN host galaxies, we rely on photo-z 
measurements. A complete description of our photo-z technique 
appears in Jouvel et al. (2013) and A. Molino et al. (2014, 
in preparation). Here, we give only a brief description of this 
technique. The spec-z and photo-z values of the SNe in our 
sample are listed in Table 2 and shown in Figure 5. 

We estimated the redshift and spectral type of all SN host 
galaxies with photometry obtained from deep Subaru images 
(in the B, V, R<,, I c , and z bands) and the Bayesian Photometric 
Redshift code (BPZ; Benitez 2000). For the host galaxies of 
SNe that were discovered in the WFC3 parallel fields, we also 
added galaxy photometry in the F125W and F160W bands. 
Some host galaxies were previously imaged by the Sloan Digital 
Sky Survey (SDSS; York et al. 2000), allowing us to include 
photometry in the u , g , r , i , and z bands as well. 

For each galaxy, BPZ calculates a likelihood, L(z, T), as 
a function of redshift, z, and spectral type, T, comparing the 


observed colors of the galaxies with the template library, and 
then multiplies it by an empirical prior, p(z, T|m), which 
depends on the galaxy magnitude in some reference band, m, 
yielding a full probability, p(z, T), for each galaxy. The new 
version of BPZ (N. Benitez 2014, in preparation) includes a 
new template library comprising six spectral energy distribution 
templates originally from PEGASE (Fioc & Rocca-Volmerange 
1997) and four early-type templates from Polletta et al. (2007). 
The PEGASE templates were recalibrated using the Wuyts et al. 
(2008) FIREWORKS photometry and spectroscopic redshifts 
to optimize its performance together with the new early-type 
galaxy templates. In total, we use five templates for early-type 
galaxies, two for intermediate galaxies, and four for starburst 
galaxies. The prior was calibrated using the GOODS-MUSIC 
(Grazian et al. 2006), Ultra Deep Field (UDF) (Coe et al. 
2006), and COSMOS (Ilbert et al. 2009) samples. Despite 
its compactness and simplicity, this library produces results 
which are comparable or slightly better than the best available 
photo-z methods (see the method comparison in Hildebrandt 
et al. 2010), which often include template libraries that are 
many times larger. As a result of the high-quality HST imaging 
used for its calibration and using an approach similar to that 
developed by Coe et al. (2006), the representation of typical 
galaxy colors provided by this library can be used to calibrate 
ground-based photometry to an accuracy of ~ 2% (A. Molino 
et al. 2014, in preparation). 

3.4. Supernova Classification 

We classify our SNe into SNe la, SNe lb/ c, or SNe Ilby fitting 
light curves to their multi-band photometry using a Bayesian ap- 
proach first introduced by Jones et al. (2013), where it was used 
to classify the CANDELS SN UDSlOWil. The full description 
of this classification technique, named the Supernova Taxonomy 
And Redshift Determination Using SNANA Templates (STAR- 
DUST), along with a detailed examination of any systematic 
biases it might introduce, will appear in a future paper by S. A. 
Rodney et al. (2014, in preparation). Briefly, for each SN we 
calculate the probability that it is a SN la, P (la), by comparing 
the observed fluxes (in all available bands and epochs) to light- 
curve models generated using the SNANA simulation package. 
We classify a SN as a SN la ifP (la) 0.5. However, as detailed 
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Figure 5. BPZ z-PDFs of the SN host galaxies. The z-PDFs are the solid curves and the spectroscopic redshifts, where available, are marked by red crosses. The 
designation of each SN appears in the upper left comer of each panel. All z-PDFs have been normalized so that P (z)dz = 1. CLF1 lDom is a “hostless” SN (see 
Section 3.5). 

(A color version of this figure is available in the online journal.) 


below in Section 4, when deriving the SN la rates, we sum the 
P (la) values of all the SNe in our sample 

The apparent magnitudes of each SN are measured with aper- 
ture photometry on the subtraction images of each epoch using 
the 1RAF routine apphot and the same apertures described in 
Section 2.1. The zero-point magnitudes and aperture correc- 
tions for ACS filters are taken from Sirianni et al. (2005). For 
WFC3-1R and WFC3-UVIS, we use the zero-point magnitudes 
calculated for a 0.4 aperture, as of 20 12 March 6, by the Space 
Telescope Science Institute. 25 The aperture corrections for the 
WFC3 filters were calculated by measuring the photometry of 
several bright stars using different apertures and adopting the 
correction for a 0. 4 aperture. 

For the SN la simulations, we use the Guy et al. (2007) 
SALT2 26 model, with nuisance parameters for the redshift, 
stretch (xl), color (c), and time of peak brightness. The core- 
collapse (CC) SNe are generated from the SNANA library 
of 43 CC SN templates, taken from the SN samples of the 
SDSS (Frieman et al. 2008; Sako et al. 2008; D’Andrea et al. 
2010), Supernova Legacy Survey (SNLS; Astier et al. 2006), 
and Carnegie Supernova Project (Flamuy et al. 2006; Stritzinger 
et al. 2009; Morrell 2012). Each of these CC SN models also 
has parameters for redshift, host extinction (A v ), date of peak 
brightness, and peak luminosity. 


25 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 

26 http://supemovae.m2p3. fr/~ guy/salt/ 


The remainder of our technique is fundamentally similar 
to other Bayesian light-curve classifiers (e.g., Kuznetsova & 
Connolly 2007; Poznanski et al. 2007a; Rodney & Tonry 2009; 
Sako et al. 2011): we compute the likelihood that a given model 
matches the observable data, multiply it by priors for the model 
parameters, then marginalize over all models to derive the final 
posterior classification probability. 

The simulated SNe are reddened and dimmed according to 
the Cardelli et al. (1989) reddening law and one of three host- 
galaxy dust extinction models: "low,” “medium,” and "high” 
dust models. For the simulated SNe la, the “low” dust model 
is the Barbary et al. (2012) skewed Gaussian fit to the Astier 
et al. (2006) SNLS, while the “mid” and “high” dust models 
are taken from Kessler et al. (2009a) and Neill et al. (2006), 
respectively. For CC SNe, we use models composed of a half 
Gaussian centered at Ay = 0 mag and an exponential of 
the form e" (Av/Tv) . These models have three parameters: the 
standard distribution of the Gaussian, a Av ; the characteristic A v 
value, Ty; and the ratio between the Gaussian and power-law 
components at A v = 0, Ao. For the “low,” “mid,” and “high” 
dust models, the values of these parameters are Oa v = 0.15, 0.6, 
and 0.5 mag; x v = 0.5, 1.7, and 2.8; and Ao = 1,4, and 3 mag. 

The peak magnitudes of each SN subtype are chosen ac- 
cording to their observed luminosity functions (LFs), which are 
detailed in Table 3. The Li et al. (2011b) LFs, derived from a lo- 
cal sample of SNe observed by the Lick Observatory Supernova 
Search (LOSS; Leaman et al. 2011; Li et al. 2011a, 2011b), 
were not originally corrected for host-galaxy extinction. Here, 
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Figure 6. Best-fitting STARDUST light curves for 15 of the 27 SN candidates. Each panel shows the multi-band photometry and best-fitting template light curve of the 
maximum-likelihood SN type. The title of each panel gives the SN name and its prior redshift, along with a legend detailing the different filters used for photometry. 
The table to the right of the light curve details the maximum-probability SN template; the posterior values of the light-curve fit parameters, including the date of peak 
brightness, redshift, and color and shape parameters; and the reduced x 2 value of the fit, x 2 - Downtumed arrows denote 3a upper limits on the flux. 

(A color version of this figure is available in the online journal.) 


we adopt “dust-free” LFs for SNe II-P and 1I-L such that when 
applying the “medium” dust model, the resultant simulated LFs 
approximate those published by Li et al. (2011b). 

A further prior is placed on the fraction of each SN type as 
a function of redshift. The distribution of SN type fractions has 
only been measured in the local universe (Li et al. 2011b), 
and it is expected to change with increasing redshift, once 
the SN la rate starts to deviate from the star-formation rate. 
Flowever, assuming a prior on the evolution of the SN type 
fraction requires us to assume a prior on the CC SN and SN la 
rates as a function of redshift. Such a prior might bias the 
SN la rates measured in this work, and so we classify our 


SN sample twice: once using this prior, and once assuming 
that the fraction of SN types remains constant with redshift. 
While the latter assumption is probably not the case in reality, 
it expresses our lack of concrete knowledge on the subject. We 
report the classification probability, P (la), of each SN in our 
sample in Table 2 both with (P (Ia) up ) and without (P (Ia) np ) the 
SN type fraction prior. The uncertainty reported for each P (Ia) wp 
value takes into account uncertainties in both the extinction and 
SN-fraction priors, while the uncertainty of P (Ia) np reflects only 
the uncertainty in the extinction prior. For each of the SNe in our 
sample, the resultant P (la) values are consistent with each other, 
and while generally P (Ia) wp < P (Ia) np , the difference between 
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Figure 7. Figure 6, continued, showing light-curve fits for the remaining 12 SN candidates. 

(A color version of this figure is available in the online journal.) 


Table 3 

SN Luminosity Functions Used for SN Classification 


Type 

Mr 

o 

Source 

la 

- 19.37 

0.47 

Wang et al. (2006) 

lb 

- 17.90 

0.90 

Drout et al. (201 1) 

Ic 

- 18.30 

0.60 

Drout et al. (201 1) 

IcBL 

- 19.00 

1.10 

Drout et al. (201 1) 

II-P 

- 16.56 

0.80 

Li et al. (2011b) 

Il-L 

- 17.66 

0.42 

Li et al. (2011b) 

Iln 

- 18.25 

1.00 

Kiewe et al. (2012) 


Notes. The Li et al. (201 lb) LFs have been corrected for host-galaxy 
extinction, as detailed in the text. 


the two values is small and has a negligible effect on the final 
SN la rates. The resultant light-curve fits, obtained without the 
SN-fraction prior, are presented in Figures 6 and 7. 

SNe that were caught on the rise and were suspected of being 
SNe la at z > 1 were followed up with further HST imaging in 
order to follow the evolution of their light curves. In our sample, 
these include CLAlOCal, CLFllVes, CLHllTra, CLP12Get, 
CLR12Arm, and CLT12Ale. Three SNe were caught sufficiently 


early, and were bright enough, to be followed up spec- 
troscopically, either from the ground or with HST. These 
were CLIllHad, CLFllVes, and CLY13Pup, whose spectra 
were obtained with Keck+LRIS, the ACS G800L grism, and 
VLT+FORS2, respectively. Using the Supernova Identification 
code (SN1D 27 ; Blondin & Tonry 2007), we classify CLI1 lHad 
and CLY13Pup as SNe la having the spec-z measured from each 
of their host galaxies. The best-fitting SN1D templates are over- 
laid on the SN spectra in Figure 2. Owing to its high redshift, 
the ACS G800L grisms caught CLF 1 1 Ves only in the rest- frame 
range -2500^1500 A, and SN1D fails to classify the spectrum 
as belonging to any type of SN at z = 1.22. When allowed to fit 
for the SN redshift, SN1D classifies CLFllVes as either a SN 
la or SN lb at z ~ 0.95. Consequently, we do not claim to have 
spectroscopic confirmation for CLF1 IVes as a SN la, although 
in the bottom panel of Figure 2, we show that the spectrum could 
be fit with the Hsiao et al. (2007) SN la template at peak and at 
z = 1.22. 

Because declining SNe appear in all four epochs, their 
photometry cannot be measured from the subtraction images. 


27 http://people.lam.fr/blondin.stephane/software/snid/ 
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CLAlONer, CLFllDom, and CLL12Luc are well offset from 
their respective host galaxies, so we measure their photometry 
from the target images and assume any contamination by 
galactic light is minimal. CLT12Ela is in a relatively faint (in the 
observed bands) area of its host galaxy. Here, too, we assume 
that any contamination by galactic light is minimal. As can 
be seen in Figure 3, CLKllBur, CLV12Gor, and CLY13Gal 
exploded in relatively bright regions of their host galaxies, so 
contamination is a certainty. However, all three of these SNe 
show signs of either a plateau or a slow decline in their light 
curves, and are classified as SNe II, as shown in Figures 6 and 7. 

3.5. Notes on Individual Supernovae 

We make the following more detailed notes on a number of 
cases with ambiguous host identification, redshift, or classifica- 
tion. 

CLDllCla has a broad z probability density function (PDF) 
with a prominent peak at z = 2.66f° 0 ° 0 7 9 and a secondary peak 
at z = 0.24;° 0 0 0 7 4 . We have attempted to obtain spectra of the 
host galaxy of this SN using several instruments (VLT+FORS2, 
Keck+LRIS, and the ACS+G800L grism). However, the host 
galaxy appears to have an absorption spectrum with no dis- 
cernible emission lines suitable for determining its redshift. At 
z ~ 2.7, the SN would be too bright to be either a SN la or 
a CC SN. While it could still be a high-z superluminous SN 
(Gal- Yam 2012), it is more likely that this is a normal CC SN 
at z = 0.24, at which its colors agree with those of CC SNe 
simulated using SNANA. In any case, we exclude this SN from 
the SN la sample. 

CLFllDom was discovered on the decline, but actually 
peaked during the period of the survey, as evidenced by its 
photometry. It has no immediately discernible host galaxy, the 
closest lying ~4.7 away, with a photo-z of 0.7k ^(, 5 . When 
classifying this SN, we chose a wide redshift prior of 0.7 ± 0.6 
and found that the light curve was best fit as belonging to a SN 
la at z = 0 . 6 , consistent with the photo-z of the nearby galaxy. 

CLL12Luc is a declining SN with four possible host galaxies. 
Of those, three have photo-z values of 0.33, 0.59, and 1.13. One 
has a spec-z of 0.36 that is consistent with its photo-z of 0.33. 
At these redshifts, the SN would be ~45, 30, and 70 kpc away 
from the core of each galaxy, respectively. Based on the colors 
of the SN, and assuming it peaked 40-70 days before the first 
observation epoch, it could either be a CC SN at z = 0.36 or 
a SN la at z “ 0.6. However, the SN colors measured from 
the last two observation epochs are inconsistent with the SN la 
solution, so in this work we treat this SN as a CC SN. 

Like CLL12Luc, CLC12Thr has several potential host galax- 
ies, for three of which we have measured the photo-z. Two 
galaxies, a face-on spiral and an elliptical, are at z = 0.23, 
while the third is an E/ SO galaxy at z ~ 1.6. At either redshift, 
the colors of the SN are consistent with those of a CC SN. How- 
ever, the z = 0.23 spiral and elliptical galaxies are much closer 
to the SN than the z ® 1.6 galaxy: -60,30, and 155 kpc, respec- 
tively. Based on these data, we classify this SN as a CC SN at 
z = 0.23. It would be interesting to inspect the elliptical galaxy, 
which is closer to the SN, for evidence of recent star formation. 

We have discovered four SNe la at z > 1.2: CLAlOCal, 
CLFllVes, CLHllTra, and CLP12Get. Although we have not 
yet succeeded in obtaining the spectroscopic redshift of the host 
galaxy of CLAlOCal, its photo-z is based on both Subaru data 
and HST NIR photometry in the F125W and F160W bands, 
so we are confident of its redshift. CLP12Get has a secure 
redshift from a photo-z fit based on nine bands from the optical 


(Subaru) to the NIR (HST). The spectrum of its host galaxy, 
obtained with VLT+X-shooter, shows no emission or absorption 
lines, but the photometry extracted from its continuum is 
consistent with the photo-z fit. CLF 1 1 Ves, as noted above, has a 
secure, spectroscopic redshift of 1.22. All four candidates have 
P (Mvp.np > 0.95. 

There are, at most, three more SNe with broad photo-z PDFs 
that could potentially enter this redshift range: CLR12Arm, 
CLT12Ale, and CLY13Gal. However, all three SNe are clas- 
sified as CC SNe according to their light curves, with pos- 
terior redshifts of -1.1, 1.4, and 0.9, respectively. CLT12Ale 
exploded in a faint host galaxy (F160W = 26.4 ±0.1 mag). 
Based on the available Subaru photometry and additional HST 
photometry in F 125W and F160W, the photo-z of the galaxy can 
only be constrained to the very wide range 0. 5-3.0. When con- 
sidering the two-parameter space spanned by the host galaxy’s 
redshift and type, we find that the most likely redshift solutions 
are either at z “ 1 or z > 5. Although superluminous SNe 
have been observed out to redshift -4, they are exceedingly 
rare (Cooke et al. 2012). It follows that the more likely redshift 
solution for this host galaxy is - 1 . This is strengthened by the 
light-curve classification of this SN, which finds a posterior red- 
shift of - 1.4. As there is some probability that CLT12Ale is a 
SN la (P (Ia) np = 0. 13 ± 0. 10), we include it in our calculation 
of the SN la rates, as detailed in Section 4, below. 

4. THE TYPE-Ia SUPERNOVA RATE 

In this section, we use the SN la sample from Section 3, 
along with the detection efficiencies as a function of magnitude 
measured in Section 3.2 and the classification probabilities 
measured in Section 3.4, to measure the rates of SNe la as a 
function of redshift, or lookback time. So as not to bias our 
results, we use the SN classification without the assumption of 
the SN-fraction prior. 

The SNe la in our sample can be divided among three 
categories, according to when they reached maximum light: 
before, during, or after the monitored interval of time spent on 
each field. Each category will have a distinct detection efficiency 
as a function of redshift. The date of maximum light can occur up 
to 40 days before and 20 days after the duration of the survey. 
These values were chosen according to the approximate time 
when the SNe la in our sample reached their peak, relative to 
the survey times in the fields where they were discovered, based 
on preliminary light-curve fits. Accordingly, the visibility time 
of our survey is defined as the sum of the times each parallel field 
in each cluster was monitored (i.e., the time between the first and 
last epoch of that field), with the addition of 40 days before and 
20 days after the observation period, in order to account for the 
SNe la in our sample that were caught either in decline or on the 
rise. Adding more or less time to the survey-extension times will 
either raise or lower, respectively, the number of SNe included 
in the rate calculation. To within Poisson errors, the change in 
extension time should, in principle, cancel out the change in SN 
numbers, leaving the resultant SN la rate unchanged. 

We define the rate, Rj a , in a redshift bin bound by redshifts 
Z\ and Z 2 , as 

Nj(z)/pi(z) 

R Ia ( Zl < z < z 2 ) = 1 ^^dv, ’ ( 3 ) 

L J ' rv J zi (1+z) dz u 

j 

where N; is the number of SNe la (see below); r|; is that 
category’s detection efficiency at the redshift, z, of each SN; 
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tj is the visibility time, composed of the time between the first 
and last epoch of observation of a field j, plus 40 days before 
the start of the survey and 20 days after its end; Aj is the solid 
angle of the searchable area of field j, divided by 4 ti steradians; 
dV are thin volume elements behind each searchable area; and 
the (1 + z) factor converts the rates from the observer frame 
to the rest frame. Although we classify a SN as a SN la if 
P (la) 0.5, we define Nj as the sum of P (Ia) np values of all the 
SNe in each subcategory (before, during, or after the monitored 
interval). This is based on our treatment of P (la) as a measure 
of the probability of a SN being a SN la. Thus, for example, 
CLD1 lCla has a P (Ia) np = 19% probability of being a SN la, 
and is counted accordingly. This approach allows us to take into 
account the uncertainty of our classifications, especially for SNe 
with sparse data. 

To compute the detection efficiency of each SN category, t|, 
we must first convert our detection efficiency from a function 
of magnitude to a function of redshift. We do this by using 
the measured detection efficiency as a function of magnitude 
from Section 3.2 to simulate the discovery process of -25,000 
SNe la. For each SN, we use the Flsiao et al. (2007) SN la 
spectral templates to simulate a light curve. The Flsiao et al. 
(2007) spectral templates are normalized so that an unredshifted, 
“normal” (i.e., with stretch s = 1) SN la at maximum light has 
a B-band apparent magnitude of M B = 0. These templates 
are first redshifted and reddened using the Cardelli et al. 
(1989) extinction law. We perform synthetic photometry on the 
redshifted and reddened spectral templates in the survey filters 
(F850LP and F160W) and construct light curves according to 
Perlmutter et al. (1999): 


l= i +M B +p - a(s - 1), 

(4) 

= t z=0 ,s=ia(l +z), 

(5) 


where m z s= i is the apparent magnitude of a “normal” SN la 
with s = 1 at redshift z; M B is the absolute magnitude of the 
SN in the B band at maximum light; p is the distance modulus 
at redshift z; a = 1.52 ± 0. 14 (Astier et al. 2006); and t s z is the 
time axis of the light curve, stretched as a result of the stretch, 
s, and time dilation at redshift z. 

Each SN is assigned a random cluster and field (WFC3 or 
ACS), redshift, M B , stretch value, and host-galaxy extinction 
value (Ay). The redshift values are drawn from a flat distribution 
in the range 0-3. Following G1 1, the absolute B-band magnitude 
at maximum light is drawn from a Gaussian centered on 
M b = - 19.37 with a standard deviation of c>m b = 0.17 (this 
standard deviation, smaller than the one used for the SN la LF 
as it appears in Table 3, reflects the spread in SN la intrinsic 
luminosity after correcting for the luminosity-stretch relation; 
Hamuy et al. 1995, 1996; Riess et al. 1996; Phillips et al. 1999). 
The stretch values, following Sullivan et al. (2006), are drawn 
from a Gaussian centered on s = 1 with a standard deviation 
of a s = 0.25 and limited to the range 0.6 < s < 1.4. This 
distribution is wide enough to account for both subluminous 
and overluminous SNe la. The distribution of the amount of 
dust in the vicinities of SNe la is as yet poorly constrained. To 
gauge the systematic uncertainty of the SN la rate caused by 
this, as in Section 3.4 above, we follow D08 and Barbary et al. 
(2012) and use four different host-galaxy extinction models in 
our simulation: Flatano et al. (1998), Riello & Patat (2005), 
Neill et al. (2006), and Kessler et al. (2009a), which have 
average extinctions of A v = 0.3, 0.3, 0.4, and 0.2 mag, 



Figure 8. SN la host-galaxy dust extinction models used in the derivation of 
the SN la rates. We use the Neill et al. (2006) model (solid black curve) as the 
fiducial model. While the models shown here go out to Ay = 7 mag, we use 
values only to Ay = 3 mag, as nearly all the objects produced by these models 
fall in the range Ay = 0-3 mag. 

(A color version of this figure is available in the online journal.) 


respectively. To remain consistent with Gil, we choose the 
Neill et al. (2006) model as our fiducial host-galaxy extinction 
model. These models are illustrated in Figure 8. Although we 
show these models (and specifically Hatano et al. 1998 and 
Riello & Patat 2005) going out to A v = 7 mag, we use values 
only out to Ay = 3 mag, after which all models produce a 
negligible number of objects. 

After choosing when the SN reached peak, the light curve 
is sampled according to the survey cadence in that particular 
cluster and field, and subtraction magnitudes are computed by 
subtracting the flux in each search epoch from the flux of the 
reference image. Using the detection efficiency at the resultant 
magnitude in each search epoch, the SN is either discovered or 
missed. The resultant detection-efficiency curves for SNe la that 
reached maximum light before, during, or after the monitored 
interval are shown in Figure 9. 

We divide the SNe in our sample into three redshift bins: 
0 < z < 0.6, 0.6 < z < 1.2, and 1.2 < z < 1.8. In each 
redshift bin, we compute the effective redshift, z e n , as 


z eff (zi < z < z 2 ) = 


Z2 

Zl 

Z2 

Zl 


— 2 — dV 

(1+Z)1 1 (Z) UV 

— - — dV 

(1+Z)T,(Z) U 


(6) 


In each redshift bin, we take the minimal and maximal 
differences between the rate as computed with the fiducial dust 
model and with the models in Figure 8 as lower and upper 
systematic uncertainties owing to dust extinction. Since we 
express the number of SNe la in our final sample as the sum of 
all the SN P (Ia) np values, we also propagate the uncertainties 
in these values and add to the rates a systematic uncertainty 
due to our classification technique. The systematic uncertainties 
from dust extinction and classification are then summed. Gil 
also considered the systematic uncertainty due to the expected 
increase in extinction as a result of dust at high redshifts (e.g., 
Mannucci et al. 2007; Mattila et al. 2012). Flowever, Gil 
did not take into account the different extinction models used 
here. Specifically, the Hatano et al. (1998) extinction model 
adds a - +9% systematic uncertainty to the SN la rate in 
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Figure 9. SN la detection efficiency as a function of redshift in the F850LP 
(solid black) and F160W (dashed red) bands for SNe la that reached maximum 
light in the B band before (bottom), during (center), and after (top) the monitored 
interval. 

(A color version of this figure is available in the online journal.) 



Figure 10. SFH measurements and parameterizations. SFH measurement 
compilations up to 2006 from HB06 are shown in filled circles, and additional 
measurements up to 20 1 0 compiled by Horiuchi & Beacom (20 1 0), here marked 
as HB 1 0, are shown as filled squares. The most up-to-date compilation from B 1 3 
is shown as open triangles. The different parameterizations are shown as curves 
and include the Cole et al. (2001) fit to the HB06 data (dashed); the Y08 (solid) 
and O08u/ 0081 (thin/ thick dotted) power-law fits; and the B 1 3 parameterization 
(dot-dashed). All data and parameterizations have been rescaled to the Bell et al. 
(2003) “diet” Salpeter IMF. 


the 1.2 < z < 1.8 bin, similar to the ~ +10% systematic 
uncertainty G1 1 added to their SN la rate at 1.5 < z < 2.0. 

While we found no SNe at z > 1.8, Figure 9 shows that 
WFC3 is still sensitive to SNe la out to z = 2.5. Consequently, 
we add a fourth redshift bin, 1.8 < z < 2.4, and compute 


Lookback Time [Gyr] 

0 12 3 4 5 6 7 8 9 10 11 



Figure 11. SN la rates from CLASH (filled, red squares) compared to rates 
from the literature and best-fitting SN la rate evolutions derived by convolving 
a power-law DTD with different SFHs. Circles denote data from surveys with 
measurements out to z ~ 1 from Cappellaro et al. (1999), Hardin et al. (2000), 
Pain et al. (2002), Tonry et al. (2003), Blanc et al. (2004), Botticella et al. (2008), 
Horesh et al. (2008), Rodney & Tonry (2010), Li et al. (201 la), Barbary et al. 
(2012), and Melinder et al. (2012). Filled circles denote the most accurate and 
precise measurements at z < 1 and are from the SDSS Stripe 82 survey (Dilday 
et al. 2010; orange), SNLS (Perrett et al. 2012; green), and SDSS DR7 (Graur & 
Maoz 2013; purple). The GOODS rates from Dahlen et al. (2008) are shown as 
triangles and the SDF rates from Graur et al. (201 1) are shown as diamonds. The 
z > 1.5 rates from these two surveys are colored in black and blue, respectively. 
The thick curves are convolutions of several SFHs (dashed, Hopkins & Beacom 
2006; solid, Yiiksel et al. 2008; thin/ thick dotted, Oda et al. 2008; dot-dashed, 
Behroozi et al. 2013) with the best-fitting power-law DTDs. The shaded area is 
the confidence region resulting from the combined 68% statistical uncertainties 
in the values of the power-law index fit with the above SFHs. The thin dashed 
lines indicate the 68% statistical uncertainty region obtained without the new 
CLASH measurements. All vertical error bars are sums of the statistical and 
systematic uncertainties. The CLASH vertical error bars are composed of the 
systematic uncertainty, shown as black thick lines, and the statistical uncertainty, 
shown as red thin lines. The horizontal error bars delineate the CLASH redshift 
bins. The Perrett et al. (2012) and z > 1.5 Dahlen et al. (2008) SN la rates have 
been shifted by Az = +0.02 to disentangle them from other results. 

(A color version of this figure is available in the online journal.) 

a 2c upper limit to the SN la rate in that bin by taking the 
95% Poisson uncertainty in the number of SNe found in the bin 
(zero), and considering the detection efficiency of the different 
SN categories in the center of the bin at z = 2.1. 

The resultant SN la rates, including both statistical and 
systematic uncertainties, are shown in Figure 11. Table 4 
summarizes the SN la rates, with and without correction for 
host-galaxy dust extinction, and Table 5 shows the complete 
error budget of our SN la rates. Table 6 compares the rates from 
this work to previous rates from the literature. Where necessary, 
the measurements have been corrected to reflect the value of 
h = 0.7 used in this work. As Perrett et al. (2012) did not take 
into account low-stretch, SN1991bg-like SNe la, we scale up 
their SN la rates by 15% (see their Section 6). As in Gil, in 
instances where rates were originally reported in units of SNuB 
(SNe per century per 10 10 L <B ; Cappellaro et al. 1999; Hardin 
et al. 2000; Pain et al. 2002; Madgwick et al. 2003; Blanc et al. 
2004), we have converted them to volumetric rates using the 
Botticella et al. (2008) redshift-dependent luminosity density 
function, 

j B (z) = (1.03 + 1.76z) x 10 8 L , B Mpc“ 3 . (7) 
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Table 4 

SN la Numbers and Rates 


Subsample 

0.0 < z < 0.6 

0.6 < z < 1.2 

1.2 < z< 1.8 

1.8 < z < 2.4 

Total 

12 

11 

4 

0 

SN host galaxies with spec-z 

10 

7 

2 

0 

Hostless SNe 

0 

1 

0 

0 

SNe la (raw) 

2.4 

6.4 

4.1 

0 

SNe la (efficiency-corrected) 

2.4 

7.0 

8.0 

0 

SN la rate without host-galaxy extinction 3 

n ,<+0.42, +0.10 
0.32,- 0.13 

n , ? +0.21,+0.12 
U - 4Z - 0.18,- 0.05 

+0.21, +0.03 
U - Z/ -0.13,-0.06 

<0.6 

SN la rate (10 -4 yr” 1 Mpc -3 ) 

n /ik+0.42,+0.10 
U,4 °- 0.32,- 0.13 

n ,,<-+0.22, +0.13 
U -^- 0.19,- 0.06 

n +0.34, +0.05 

0.22,- 0.09 

< 1.7 

Effective redshift 

0.42 

0.94 

1.59 

2.1 


Notes. The 1.8 < z < 2.4 rate is a 2o upper limit. 

3 The errors, separated by commas, are respectively the 68% Poisson statistical uncertainties on the number of SNe, and systematic uncertainties 
due to possible misclassification and different host-galaxy extinction models, respectively. 


Table 5 

SN la Rate Uncertainty Percentages 


Uncertainty 

0 < z < 

0.6 

V 

so 

o 

?< 1.2 

1.2 < z < 

: 1.8 

Poisson 

+92/ - 

71 

+49/ 

- 44 

+75/ - 

49 

Classification 

+ 17/ - 

28 

+4.6/ 

- 9.4 

+ 1.9/ - 

2.7 

Dust model 

+5.4/ - 

0.1 

+25/ 

- 3.3 

+8.6/ - 

18 

Total 

+114.4/ - 

99.1 

+78.6/ 

- 56.7 

+85.5/ - 

69.7 


Note. All uncertainties are reported as percentage of the rates. 


5. THE TYPE-Ia SUPERNOVA DELAY-TIME 
DISTRIBUTION 

In this section, we test different models of the DTD by 
convolving them with various cosmic SFHs and fitting the 
resultant SN la rate histories to the SN la rate measurements 
from the previous section, along with rates from the literature. 
We include all the rate measurements from Table 6 except 
for Neill et al. (2006, 2007), which have been superseded by 
Perrett et al. (2012); Dahlen et al. (2004) and Kuznetsova et al. 
(2008), which have been superseded by D08; Barris & Tonry 
(2006), which has been superseded by Rodney & Tonry (2010); 
Poznanski et al. (2007b), which has been superseded by Gil; 
and Madgwick et al. (2003) and Krughoff et al. (2011), which 
have been superseded by Graur & Maoz (2013). We do not use 
thez > 2 upper limit from the previous section as it is too high to 
affect the DTD fits. In total, we use 50 SN la rate measurements, 
of which 4 1 are at z < 1 and 9 are at z > 1 . 

Following Gil, we test different SFHs, including the Cole 
et al. (200 1) parameterization fit to the data collected by Hopkins 
& Beacom (2006, HB06); the SFH presented by Yiiksel et al. 
(2008, Y08) and upper (O08u) and lower (0081) limits from Oda 
et al. (2008, 008) which can be approximated as broken power 
laws with a break at z = 1 and with varying indices before and 
after the break; and the recent Behroozi et al. (2013, B13) SFH. 
These SFHs, and the data they are based on, are presented in 
Figure 10. 

When deriving SFH measurements, various authors use 
different versions of the initial-mass function (IMF), leading to 
different scalings of the SFH. In order to maintain consistency 
across the different SFHs, we must choose one IMF and rescale 
the SFHs accordingly. As in G1 1 and Graur & Maoz (2013), we 
assume a “diet” Salpeter IMF (Bell et al. 2003), which is similar 
to the Salpeter (1955) IMF with lower and upper mass limits of 
0.1 and 125 M , respectively, but with a stellar mass-to-light 
(M/L) ratio that is scaled down by a factor of 0.7 in order to 


fit the M/L ratios measured in disks (Bell & de Jong 2001). 
The choice of this IMF requires us to scale down the SFHs of 
HB06, 008, and Y08, who assumed a Salpeter (1955) IMF, by a 
factor of 0.7. The B13 SFH, where a Chabrier (2003) IMF was 
assumed, is scaled up by a factor of 0.7. To allow comparisons 
with Gil, we use the Y08 SFH as the fiducial model in our 
DTD recoveries and the other SFHs to estimate a systematic 
uncertainty in the values of the parameters of the DTD model 
tested below. 

We test a power-law DTD of the form 'Pit) = 
^iGyr(t/ 1 Gyr)P, setting its index, p, and scaling, ^FiGyo as 
free parameters, leaving 48 degrees of freedom for the fit. 
The DTD is set to zero before 40 Myr, to allow for 8 M 
stars to evolve into CO WDs. The Y08 SFH yields a best- 
fit index value of P = - l.OO+o o^m) with a reduced x 2 
(Xv) 0-7, where the statistical uncertainties are the 68% 
and 95% (in parentheses) confidence regions, respectively. The 
other SFHs yield a systematic uncertainty of tg'og, with x 2 
values in the range 0.7-0. 8, yielding a final value of P = 
- J-OO-op^oTo) (statistical)^ 1 ^ (systematic). This value is con- 
sistent with those obtained by Gil and Graur & Maoz (2013) 
and in a variety of different SN surveys and using different 
DTD recovery techniques (see Maoz & Mannucci 2012). Inte- 
grating the DTD over a Hubble time, we find that the num- 
ber of SNe la per formed mass, N/M*, lies in the range 
(0.5-1. 3) x 10~ 3 SNe M -1 , similar to the ranges found by 
Gil and Graur & Maoz (2013). The best-fitting SN la rate his- 
tories derived from each SFH, along with the 68% uncertainty 
region, are shown in Figure 11. 

We investigate also the viability of a Gaussian DTD fit 
(Strolger et al. 2004, 2010) to the SN la rates. We start by 
testing the Gaussian DTD proposed by D08, with a mean delay 
time of 3.4 Gyr and a standard deviation of 0.68 Gyr. As in Gil, 
we allow the scaling of this DTD to vary as a free parameter. 
Coupled with the SN la rates, the only SFH that does not disfavor 
this DTD is the HB06 SFH, with a reduced x 2 value of 1.1. All 
other SFHs result in SN la rate evolutions that are excluded 
at a >95% significance level, with the 008 and B13 SFHs 
specifically excluded at a > 99% significance level. We next test 
a general Gaussian DTD, where we allow the mean delay time, 
standard deviation, and scaling to vary as free parameters, while 
requiring that 95% of the area under the DTD remain above a 
delay time of 40 Myr (thus ensuring the resultant DTD retains 
a Gaussian shape). The B13 and lower-limit 008 SFHs result 
in Gaussians that are excluded at a > 95% significance level. 
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Table 6 

SN la Rate Measurements 


Redshift 

Nia 

Rate (10 4 yr 1 Mpc 3 ) 

Reference 

Redshift 

Nia 

Rate (10 4 yr 1 Mpc 3 ) 

Reference 

0.01 

70 

0.183 ± 0.046 

Cappellaro et al. (1999) a 

0.55 

72 

n ^+0.07, +0.05 
U --°- 0.07,- 0.06 

Perrett et al. (2012) e 

<0.019 

274 

rv r\s ^+0.034, +0.043 
U.ZOj_ 0 Q 33 _ 0 043 

Li et al. (201 la) b 

0.552 

41 

/-v £0+0. 10, +0.26 

U - OJ -0.10,-0.27 

Neill et al. (2007) d 

0.0375 

516 

O 970+0-112, +0.015 
U,Z/5 - 0.083,- 0.000 

Dilday et al. (2010) c 

0.62 

7 

, ?q +0.88,+0.27 
lZy - 0.57,- 0.28 

Melinder et al. (2012) 

0.09 

17 

0 dq+ 0.09 
uzy - 0.07 

Dilday et al. (2008) 

0.65 

23 

1.49 ± 0.31 

Barris & Tonry (2006) d 

0.098 

19 

0 24 +012 
u - z ^- 0.12 

Madgwick et al. (2003) a,d 

0.65 

10.09 

0 4Q+0-17,+0.14 
u -^ y -0.17,-0.08 

Rodney & Tonry (2010) 

0.1 

516 

7cq+0-052,+0.018 
U - ZDV - 0.044,- 0.001 

Dilday et al. (20 1 0) c 

0.65 

91 

n <-^+0.06, +0.05 
0.06,- 0.07 

Perrett et al. (20 12) e 

0.1 

52 

O caq+0.098, +0.058 
U - 3 -0.085,-0.047 

Rrughoff et al. (20 1 1 ) d 

0.714 

42 

, ,o+0. 19, +0.54 
A,A ' 5 - 0.19,- 0.70 

Neill et al. (2007) d 

0.11 

90 

O 747+0.029,+0.016 
U - Z4 '-0.026,-0.031 

Graur & Maoz (2013) 

0.74 

5.5 

Q 4 +0.36 
0.32 

Poznanski et al. (2007b) d 

0.13 

14 

n , ™+0.056, +0.035 
y . 135_ 0043 ,- 0.035 

Blanc et al. (2004) a 

0.74 

20.3 

A 70+0.33 
u - y - 0.41 

Graur et al. (2011) 

0.14 

4 

n ^o+0.22, +0.07 
U - ZS - 0.13,- 0.04 

Hardin et al. (2000) a 

0.75 

28 

1.78 ± 0.34 

Barris & Tonry (2006) d 

0.15 

516 

n 9n7 +0 038 ’ +0 - 035 
U.3U /_ o.034,- 0.005 

Dilday et al. (2010) c 

0.75 

14.29 

n £g+0.21, +0.23 
11.05-021,-0.14 

Rodney & Tonry (2010) 

0.15 

1.95 

77 +0.23,+0.07 
U - JZ - 0.23,- 0.06 

Rodney & Tonry (2010) 

0.75 

110 

/-, £-7+0. 07, +0.06 
U - 0/ - 0.07,- 0.08 

Perrett et al. (2012) e 

0.16 

4 

/-v i £+0. 10, +0.07 
U1 °- 0.10,-0.14 

Perrett et al. (2012) e 

0.80 

14 

, S7 +0.44,+0.75 
L - D/ - 0.25,- 0.53 

Dahlen et al. (2004) d 

0.2 

17 

0 189 +0 042,+0 ' 018 +0 42 
u - loy - 0.034,- 0.015 z u -^ z 

Horesh et al. (2008) 

0.80 

18.33 

a 90+0.25 
uv:) -0.25 

Kuznetsova et al. (2008) d 

0.2 

516 

n 7/10+0.032, +0.082 
U.345_ 0.030,- 0.007 

Dilday et al. (20 1 0) c 

0.807 

5.25 

, , 0+O.6O, +0.44 
A,1 °- 0.45,-0.28 

Barbary et al. (2012) 

0.25 

1 

0.17 ± 0.17 

Barris & Tonry (2006) d 

0.83 

25 

1 0 a+0.33,+0.73 
i.JU_ 0 27,-o.51 

Dahlen et al. (2008) 

0.25 

516 

O 7£c+0.031,+0.182 

\J.jOD_ 0.028,- 0.012 

Dilday et al. (2010) c 

0.85 

15.43 

n 7 2+0.22,+0.31 
U - 0.22,- 0.16 

Rodney & Tonry (2010) 

0.26 

16 

n t^+0.08,+0.07 
U " 3Z - 0.08,-0.08 

Perrett et al. (2012) e 

0.85 

128 

n ££+0.06,+0.07 
U - 00 - 0.06,- 0.08 

Perrett et al. (20 12) e 

0.3 

31.05 

-5^j+0. 16, +0.21 

U - J4 - 0.15,-0.22 

Botticella et al. (2008) f 

0.94 

6.4 

Q , 5 +0.22,+0.13 
u -^ J -0.19,-0.06 

CLASH (this work) 

0.3 

516 

rv /i 74+0-037, +0.396 
U.+3+- 0.034,-0.016 

Dilday et al. (2010) c 

0.95 

13.21 

n 7 £+0.25,+0.32 
U - /0 - 0.25,- 0.26 

Rodney & Tonry (2010) 

0.35 

5 

0.530 ± 0.024 

Barris & Tonry (2006) d 

0.95 

141 

n on+0.09,+0. 12 
y, oy - 0.09,- 0.14 

Perrett et al. (2012) e 

0.35 

4.01 

n -,,+0.19, +0.07 
U - J4 - 0.19,-0.03 

Rodney & Tonry (2010) 

1.05 

11.01 

n 7 q0.28,+0.36 
U - /y - 0.28,- 0.41 

Rodney & Tonry (2010) 

0.35 

31 

,, +0.07, +0.06 
U -^ A - 0.07,- 0.07 

Perrett et al. (2012) e 

1.05 

50 

n 8 c+0.14,+0.12 

U.OD-o 14-0.15 

Perrett et al. (2012) e 

0.368 

17 

n , +0.05,+0.08 
0.05,- 0.03 

Neill et al. (2007) d 

1.187 

5.63 

1 0 9+0. 65, +0.69 

A -DD- O.49, -0.26 

Barbary et al. (2012) 

0.40 

3 

n .q+0.34, +1.54 
0.27,-0.25 

Dahlen et al. (2004) d 

1.20 

6 

, , ^+0.47, +0.32 
0.26,- 0.44 

Dahlen et al. (2004) d 

0.40 

5.44 

0 53 +0 ' 39 

U.DD-o 17 

Kuznetsova et al. (2008) d 

1.20 

8.87 

0 75 +0 - 35 
u - /:> -0.30 

Kuznetsova et al. (2008) d 

0.42 

2.4 

n ,£+0.42, +0.10 
°- 46 - 0.32,-0.13 

CLASH (this work) 

1.21 

20 

, 77 +0.36,+0.38 
A - JZ - 0.29,- 0.32 

Dahlen et al. (2008) 

0.442 

0 

n nn +0.50,+0.00 
UUU - 0.00,-0.00 

Barbary et al. (2012) 

1.23 

10.0 

1-05^ 

Poznanski et al. (2007b) d 

0.45 

9 

0.73 ± 0.24 

Barris & Tonry (2006) d 

1.23 

27.0 

0-84!° d f 8 

Graur et al. (201 1) 

0.45 

5.11 

n 71 +0.15, +0.12 
U - M - 0.15,- 0.04 

Rodney & Tonry (2010) 

1.535 

1.12 

n 77 + 1.07, +0.44 
U - ' '-0.54,-0.77 

Barbary et al. (2012) 

0.45 

42 

n ,,+0.07, +0.05 
U '^ A - 0.07,- 0.06 

Perrett et al. (2012) e 

1.55 

0.35 

n 12 +0 - 58 
U - 1 Z - 0.12 

Kuznetsova et al. (2008) d 

0.46 

8 

0.48 ± 0.17 

Tonry et al. (2003) 

1.59 

4.1 

n /i^+oj^+o.os 
0.22,- 0.09 

CLASH (this work) 

0.467 

73 

, 7 +0.06,+0. 13 
U -^ Z - 0.06,- 0.09 

Neill et al. (2006) d 

1.60 

2 

Q 44 +0.32, +0.14 
- 0.25,-0.11 

Dahlen et al. (2004) d 

0.47 

8 

H on+0-37,+1.66 
U.5U- 0.27,-0.26 

Dahlen et al. (2008) 

1.61 

3 

n 47 +0.39,+0.19 
u -^ z - 0.23,- 0.14 

Dahlen et al. (2008) 

0.55 

38 

H c/ro+0.098, +0.098 

U.305- o. 0 88,- 0.088 

Pain et al. (2002) a 

1.67 

3.0 

a 01 +0.79 
U - 81 -0.60 

Poznanski et al. (2007b) d 

0.55 

29 

2.04 ± 0.38 

Barris & Tonry (2006) d 

1.69 

10.0 

1 02 +0 - 54 
a - uz -0.37 

Graur et al. (2011) 

0.55 

6.49 

r\ 77 +0. 14, +0.07 
U - JZ - 0.14,-0.07 

Rodney & Tonry (2010) 

2.1 

0 

< 1.7 

CLASH (this work ) 8 


Notes. Redshifts are means over the redshift intervals probed by each survey. Ni a is the number of SNe la used to derive the rate. Where necessary, rates have been 
converted to h = 0.7. Where reported, the statistical errors are followed by systematic errors, and separated by commas. Rates from this work are shown in bold. 
a Rates have been converted to volumetric rates using Equation (7). 
b Li et al. (2011a) consider SNe la within 80 Mpc. 
c Dilday et al. (2010) compute their rates using 516 SNe la at z < 0.5. 
d These measurements have been superseded by more recent results, as detailed in Section 5. 

e Perrett et al. (2012) do not include SN1991bg-like SNe la in their rates. Here, their measurements are scaled up by 15% (see their Section 6 ). 
f Botticella et al. (2008) found a total of 86 SN candidates of all types. See their Section 5.2 for details on their various subsamples and classification techniques. 

8 2a upper limit on the SN la rate, as derived in Section 4. 


The Y08, upper-limit 008, and HB06 SFHs, on the other hand, 
result in Gaussians with means in the range p = 2. 7-3. 3 Gyr 
with standard deviations of a = 0.8-1. 6 Gyr and reduced y_ 2 
values of 0.9-1. 3. These Gaussian DTDs are centered at slightly 


lower mean delay times, but are wider, than the D08 Gaussian 
DTD. The resultant fits to the SN la rate evolution are shown 
in Figure 12. Although at first sight, it might appear that the 
z > 1.5 SDF rate measurement is driving the exclusion of the 
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Red shift 


0.0 0.5 1.0 1.5 



Red shift 



Figure 12. Observed SN la rates compared to predicted SN la rate evolution 
tracks from the convolution of different SFHs with a best-fitting (top) Gaussian 
DTD and (bottom) broken power-law DTD of the form ^(t) = H^t/ 1 Gyr) - 1/2 
up to t c = 1.7 Gyr, and 4^(t) oc t^ afterward. Symbols are as in Figure 11. 
The SFHs used for each fit are listed in the top panel along with best-fitting 
parameter values: the mean, p , and standard deviation, a , of the Gaussian DTD, 
in Gyr; and the slope of the second power law, P , of the broken power-law DTD. 
The CLASH upper limit at z > 1.8, which was not used in the fits, has been 
removed. 

(A color version of this figure is available in the online journal.) 


Gaussian DTDs, most of the fitting power actually comes from 
the accurate z < 1 SDSS and SNLS measurements. 

As in G1 1, we test the further possibility that at early times the 
DTD is dominated by the production efficiency of double WD 
systems, which is described by a power law of the form t -0 5 
(Pritchet et al. 2008), until some cutoff time, t c , when a second 
physical process takes over, described as a power law having 
a different slope, t^. Whereas in Gil, we set the slope of the 
second power law to (3 = - 1 and fit for t c , here we set t c = 1.7 
Gyr, the lifetime of a 2 M star, the least massive star expected 
to produce a ~0.7 M CO WD (see Figure 4 of Girardi et al. 
2000), and find that the best-fitting slope of the second power 
law is P = - 1 ■46_° 0 l | f 3 ( C J ) 2 2 6 2 ) ) (statistical) _ 0 0 48 (systematic) with 
reduced % 2 values in the range 0.7-0. 9. 


Finally, we test both DD and SD DTDs resulting from binary 
population synthesis (BPS) simulations. Flere, we use updated 
versions of the scaled models presented in Figures 2 and 3 
of Nelemans et al. (2013; for color versions of these DTD 
figures see Wang & Flan 2012; updated versions of the models 
courtesy of G. Nelemans 2013, private communication). As in 
Nelemans et al. (2013), we designate the BPS DTD models 
by the groups that computed them: Yungelson (e.g., Yungelson 
2010), the Yunnan group (Wang/ Flan et al.; e.g., Wang et al. 
2010), the StarTrack code (Ruiter et al.; e.g., Ruiter et al. 2009), 
the Brussels group (Mennekens et al.; e.g., Mennekens et al. 
2010), the Utrecht group (Claeys et al.; e.g., Claeys et al. 2013), 
and the SeBa code (Bours/ Toonen; e.g., Toonen et al. 2012). 
We present the updated versions of the scaled BPS DTD models 
in Figure 13. As Nelemans et al. (2013) scaled the different 
DTDs to a Kroupa et al. (1993) IMF, we rescale them to the 
diet Salpeter IMF by multiplying them by a factor of 0.7. For 
comparison with the volumetric SN la rates, we convolve the 
different BPS DTDs with the B13 SFH, as parameterized by 
their Equation (FI), which we reproduce here as 

CSFR(z) = — — — — -, (8) 

V ' IQA(z-zo) + JQB(z-zo) v ' 

where CSFR(z) is the cosmic SFH as a function of redshift, and 
the constants A B, C, and z 0 are given in Table 7 of B13 as 
A = - 0.997, B = 0.241, C = 0.180, and z 0 = 1.243. 

As has been commented elsewhere (e.g., Maoz & Mannucci 
20 1 2; N elemans et al. 20 1 3), the BPS DTDs, both for DD and SD 
scenarios, fail to produce the number of observed SNe la. Here, 
we test the shape of the BPS DTDs by treating their scaling 
as a free parameter. The resultant SN la rate evolutions are 
presented in Figure 14. The BPS DD models require scalings by 
small factors of 3-9 and result in reduced x 2 values of/, 2 1, 

consistent with the SN la rates. The BPS SD models, on the 
other hand, require large scaling factors of > 10 (except for 
the Wang/ Han and Mennekens DTDs, which require scaling 
factors of -4) and result in reduced /_ 2 values ofy_ 2 > 1.8, thus 
excluding all BPS SD models at a > 99% significance level. 

The poor fits of the SD models are the result of their low 
DTD amplitudes at long delay times. The accurate and precise 
z < 1 SN la rates, which are most sensitive to the long 
delay-time component of the DTD, have the most leverage on 
the scalings of the BPS SD DTDs. Because the SD models 
have low amplitudes at long delay times, the z < 1 SN la rates 
force scalings of large factors that then cause the resultant SN 
la rate evolutions to overshoot the z > 1 rates. 

It is instructive to compare some of the SD DTD models, 
and their resultant SN la rate evolutions, in detail. The Claeys 
model has the lowest amplitude at long delay times, which is 
why its resultant SN la rate evolution overshoots all the other 
models in the bottom panel of Figure 14. On the other hand, 
although the Yungelson model only has an intermediate delay- 
time component, that component is at longer delay times than 
the Claeys model, so it results in a SN la rate evolution with 
less amplitude than the rate evolution produced by the Claeys 
model. Finally, the Ruiter model has the highest amplitude at 
long delay times, after the Wang/ Han model, which results in a 
low scaling. However, the Ruiter model has lower amplitude 
at short delay times, compared to all other DTD models 
besides Yungelson, which is why the low scaling, forced by 
the long-delay component, results in a SN la rate evolution that 
undershoots the z > 1 SN la rates. 
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Delay time [Gyr] 


Figure 13. Binary population synthesis DTD models for the DD (top) and SD 
(bottom) scenarios, compared to observations. The BPS DTD models, shown 
here as different-colored curves, are updated versions of the ones that appear 
in Nelemans et al. (2013; G. Nelemans, private communication) and are listed 
according to the groups that computed them: Yungelson (e.g., Yungelson 2010; 
black solid), the Yunnan group (Wang/ Han et al.; e.g., Wang et al. 2010; 
blue dashed), the StarTrack code (Ruiter et al.; e.g., Ruiter et al. 2009; red 
dot-dashed), the Brussels group (Mennekens et al.; e.g., Mennekens et al. 2010; 
green solid), the Utrecht group (Claeys et al.; e.g., Claeys et al. 2013; orange 
dot-dashed), and the SeBa code (Bours/Toonen; e.g., Toonen et al. 2012; purple 
dashed). For comparison, we also show reconstructed components of the DTD 
from observations of SNe la in 0.4 < z < 1.2 elliptical galaxies (Totani 
et al. 2008; white diamonds), galaxy clusters (Maoz et al. 2010; blue triangles), 
LOSS-SDSS galaxies (Maoz et al. 2011; green circles), SDSS-II galaxies (Maoz 
et al. 2012; red right-pointing triangles), and SDSS DR7 spectra (Graur & Maoz 
2013; orange squares). 

(A color version of this figure is available in the online journal.) 


6. CONCLUSIONS 

In this work, we have presented a sample of 27 SNe discov- 

ered in the parallel fields of the 25 CLASH galaxy clusters. 
Of these, ~ 13 were classified as SNe la, four of which are at 
z > 1.2. Using the SN la sample, we measured the SN la rate out 

to z ~ 1.8 and obtained an upper limit on the rate in the redshift 
range 1.8 < z < 2.4. Within the uncertainties of all the mea- 
surements, these rates are consistent with both the HST/ GOODS 
and the Subaru/ SDF SN la rates. Based on these rates, along 
with previous rates from the literature, we have shown that when 
convolved with different cosmic SFHs, a power-law DTD with 


Redshift 


0.0 0.5 1.0 1.5 




Figure 14. Observed SN la rates compared to predicted SN la rate evolution 
tracks from the convolution of the Behroozi et al. (2013) SFHs with BPS DD 
(top) and SD (bottom) DTD models from the literature (Nelemans et al. 2013). 
Symbols are as in Figure 12, except for the SN la rate evolution curves, which 
are labeled in the top panel according to the DTD used in each fit, as shown in 
Figure 13. 

(A color version of this figure is available in the online journal.) 


an index of- LOO/° 0 ° 0 6 * * /!J ) () 1 9 o) ( statistical)^ * 0 g (systematic) is con- 
sistent with the data. The systematic uncertainty derives from 
the wide range of possible SFHs considered. 

We have also shown that the overall shape of DTDs from BPS 
DD models are consistent with the SN la rate measurements, 
as long as the models are scaled up by factors of 3-9, while all 
BPS SD models are ruled out at a > 99% significance level. 

The SN la rates at z < 1 require a DTD with a significant 
delayed component, such as the power-law DTD tested here. 
The high-redshift SN la rates provide a probe of the early times 
of the DTD, where the DTD could either continue with an index 
of ~ - 1, as found here, or perhaps transit to a lower index of 
- 0.5, as proposed by Pritchet et al. (2008). If the SD scenario 
contributes significantly to the SN la rate, as claimed by some 
recent work (e.g., Sternberg et al. 2011 and Dilday et al. 2012), 
its main effect would be on the high-redshift rates. However, 
to have any real discriminatory power on the different DTD 


17 


The Astrophysical Journal , 783:28 (19pp), 2014 March 1 


Graur et al . 


models, the SN la rates at z > 1 must be more accurate and 
precise than they currently are. To make the most efficient use 
of the CLASH SN sample, we will combine it with the final 
CANDELS sample in a future paper. Together, the two samples 
will contain a similar number of SNe la as the Subaru/ SDF 
sample from Gil. However, their systematic uncertainties will 
be lower, as they will make use of light curves and spectroscopy, 
where available, as done in this work. Finally, the upcoming HST 
Frontier Fields program- ' (PI: M. Mountain) will observe six 
pairs of galaxy clusters and blank fields containing field galaxies, 
using 140 orbits of ACS and WFC3 for each pair of galaxy 
cluster/ blank field during Cycles 21-23. Based on our work 
on CLASH and CANDELS, we expect that this survey, which 
will go deeper than either of the previous surveys, will discover 
~20 SNe, including five z > 1.5 SNe la. Once this sample 
is added to the combined CLASH+CANDELS SN sample, we 
may finally have high-redshift SN la rates accurate enough to 
probe the early part of the DTD. 
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